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^ ' Viable alternatives to astrophysical black holes include hyper-compact objects without horizon, 

^ . such as gravastars, boson stars, wormholes and superspinars. The authors have recently shown that 

. typical rapidly-spinning gravastars and boson stars develop a strong instability. That analysis is 

extended in this paper to a wide class of horizonless objects with approximate Kerr-like geometry. 
A detailed investigation of wormholes and superspinars is presented, using plausible models and 
mirror boundary conditions at the surface. Like gravastars and boson stars, these objects are 
unstable with very short instability timescales. This result strengthens previous conclusions that 
observed hyper-compact astrophysical objects with large rotation are likely to be black holes. 
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PACS numbers: 04.40.Dg,04.30.Nk,04.25.Nx,95.85.Sz,04.70.-s 



I. INTRODUCTION 



Astrophysical Black Holes (BHs) arc believed to be common objects in galaxies. Their mass is expected to span 
. many orders of magnitude, from a fraction of the solar mass (primordial BHs in the galactic halo) to few solar masses 
' (stellar BHs in the galactic plane) up to several billions of solar masses (supermassive BHs in galactic centers) . Their 
t-H , angular momentum should be close to the extremal limit due to accretion and mergers [l], Q • For example, if quasars 
are powered by supermassive BHs, astrophysical observations suggest that they should be rotating near the Kerr 
^ ! bound [|. 

Unquestionable observational evidence of the existence of BHs is still lacking [1, d, Q . Current astrophysical data 
cannot rule out "BH forgeries" , i.e. hyper-compact objects with redshift and geodesies similar to those of BHs, but 
lacking an event horizon. Several models of hyper-compact objects with these characteristics have been known in 
the literature for some time. Among these models, gravastars 0, Q and boson stars [H, [l(| have been proposed as 
the most viable alternatives to astrophysical BHs. The authors recently showed that rapidly spinning gravastars and 
boson stars may develop a strong ergoregion instability [Tlj |. Their typical instability timescales are of order of 0.1 
seconds to 1 week for objects with mass M = 1 — 10 6 A/ Q . Therefore, observed astrophysical hyper-compact objects 
■ are likely not to be gravastars nor boson stars. 

The purpose of this paper is to compute the ergoregion instability for other horizonless, Kerr-like hyp er-compact 



objects: wormholes and superspinars. Wormholes can be objects even simpler than BHs [12[ Il3l 1 1-411 - They are 



infinitesimal variations of the Schwarzschild space-time which may be indistinguishable from BHs [151 ]. In a string 
theory context, the fuzzball model replaces BHs by horizonless structures [16] . The BH-like geometry emerges in 
a coarse-grained description which "averages" over horizonless geometries and produces an effective horizon at a 
radius where the individual microstatc geometries start to differ. Superspinars are solutions of the gravitational field 
equations that violate the Kerr bound. These geometries could be created by high energy corrections to Einstein 
gravity such as those present in string- inspired models (TtI [Lsj . Superspinars are expected to have compactness of the 
order of extremal rotating Kerr BHs and to exist in any mass range. 
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A rigorous analysis of the ergoregion instability for these models is a non-trivial task; known wormhole solutions 
are special non- vacuum solutions of the gravitational field equations. Thus their investigation requires a case-by-case 
analysis of the stress-energy tensor. Exact solutions of four-dimensional superspinars are not known. To overcome 
these difficulties, the following analysis will focus on a simple model which captures the essential features of most 
Kerr-likc horizonlcss hyper-compact objects. Superspinars and rotating wormholes will be modeled by the exterior 
Kerr metric down to their surface, where Dirichlet boundary conditions are imposed. This problem is very similar 
to Press and Teukolsky's "BH bomb" [SEl, i.e. a rotating BH surrounded by a perfectly reflecting mirror with 
its horizon replaced by a reflecting surface. These boundary conditions are perfect mirror conditions and require a 
reflection coefficient R = 1 . In a more realistic model R < 1 and a certain transmittance T = 1 — R should be taken 
into account, which will in principle decrease the strength of the ergoregion instability. We argue that the qualitative 
behavior of the instability is the same as long as the reflection or the superradiant amplification are large enough. 
Letting p be a superradiant factor, one expects an ergoregion instability to develop whenever p(l — T) > 1. In the 
perfect mirror limit T = and the superradiant condition is simply p > 1. The general case can be handled using 
both the analytical and numerical techniques presented here. 

In Section [ill we introduce the class of objects we will deal with in this work. They are general approximations 
to superspinars and wormholes with the basic key features retained. In Section IIIII we show how to solve for the 
instability analytically in two different regimes. The details of these computations are left for Appendices [A] and [B] 
These approximations are compared with numerical results in Section IIV1 where we also show that another kind of 
instability sets in for general naked singularities. This "algebraic" instability can be computed algebraically in the 
Kerr geometry. We close with a brief discussion of our results. 

II. SUPERSPINARS AND KERR-LIKE WORMHOLES 

A superspinar of mass M and angular momentum J = aM can be modeled by the Kerr geometry [l7| 



<l»i, ,., = - ( 1 - KA dt 2 + ^-dr 2 



'Kerr 

where 



A 



(r 2 + a 2 ) sin 2 9 + ^-^a 2 sin 4 i 



4Mr 

asin 2 9dcj)dt + me 2 , (2.1) 



E = r 2 + a 2 cos 2 9, A = r 2 +a 2 -2Mr. (2.2) 



Unlike Kerr BHs, superspinars have a > M and no horizon. Since the domain of interest is — oo < r < +oo, the 
space-time posesses naked singularities and closed timelike curves in regions where < (see the monograp h by 
Chandrasekhar [34[ ). High energy modifications in the vicinity of the singularity are also expected. Following Ref. [l7j . 
a small region around the origin is excised or assumed to be modified by, say, stringy corrections. The most popular 
excision method uses domain walls formed by supertubes fl7l [2lj . Kerr-like wormholes are described by metrics of 
the form 

<^ s wormholc = ^ s Kcrr + Sg a bdx a dx b , (2-3) 

where 5g a b is infinitesimal. In general, Eq. I|2.3p describes an horizonless object with excision at some small distance 
of order e from the would-be horizon. (See Ref. [l5| for details on nonrotating wormholes). A detailed consideration 
of rotating stationary wormholes, throat location and conditions on the metric is given by Teo [12] to which we refer 
for further details. Here, we simply assume these conditions are satisfied. Wormholes require exotic matter and/or 
divergent stress tensors, thus some ultra-stiff matter is assumed close to the would-be horizon. In the following, both 
superspinars and wormholes will be modeled by the Kerr metric with a rigid "wall" at finite Boyer-Lindquist radius 
r , which excludes the pathological region. We will consider both a/M < 1 and a/M > 1. 

If the background geometries of superspinars and wormholes are sufficiently close to the Kerr geometry, their 
dynamical perturbations are determined by the equations of perturbed Kerr BHs. The proof is straightforward. 
Consider a minimally-coupled scalar field <I> propagating on a space-time with metric 

9ab = g° ab + 6g a b , (2.4) 



where 6g a b <C g a b- At first order in 8g a b, the Klein-Gordon equation reads 



d a 



-</;/"'wi> ( i + \g cd 5gc d 



d a 



-</<V/""<-Vl> ! 1 + \g cd 5g c 



= 0. (2.5) 



If g cd 5g c d <§C 1, Eq. (|2.5p is identical to the equation of a scalar field in the Kerr geometry. This result also generalizes 
to the Maxwell field. Gravitational perturbations can be handled as in the Kerr geometry only if they are larger than 
5g a b at any time. 
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III. INSTABILITY ANALYSIS: ANALYTIC RESULTS 



The instability of superspinars and wormholcs is studied by considering Kerr geometries with arbitrary rotation 
parameter a and a "mirror" at some Boyer-Lindquist radius r . Using the Kinnersley tetrad and Boyer-Lindquist 
coordinates, it is possible to separate the angular variables from the radial ones, decoupling all quantities. Small 
perturbations of a spin-s field are reduced to the radial and angular master equations [23| 

A~ s — { A s+1 dRlm 
dr \ dr 



K z - 2is(r - M)K 
A 



+ Aiscur — A 



o. 



(3.1) 



[(1 - x 2 ) s Si mtX ] + 



(acox) 2 — 2au>sx + s + s Ai Tl 



(m + sx) z 



(3.2) 



where x = cos9, A = r 2 — 2Mr + a 2 , K = (r 2 + a 2 )uj — am and the separation constants A and s Mm are related by 



A = s A lm + a uj 



2, ,2 



2amLU . 



(3.3) 



If a ^ M, the space-time possesses one or two horizons located at r± = M ± V 'M 2 — a 2 . Equations (|3.ip - (|3.2p can be 
analytically solved in the slowly-rotating and low- frequency regime, u)M <C 1 [UJ [H, [2j| , and in the rapidly-spinning 
regime, where r + ~ r_ and ui ~ nin^, = a/(2Mr+) being the angular velocity at the horizon (25j . 



A. Slowly rotating objects and low frequencies 




FIG. 1: Imaginary and real parts of the characteristic gravitational frequencies for an object with a — 0.998M, according to 
the analytic calculation for rapidly-spinning objects. The mirror location is at ro = (1 + e)r+. The real part is approximately 
constant and close to mfl, in agreement with the assumptions used in the analytic approach. 



The instability timescales for small rotation and low frequencies can be computed by approximating the Teukolsky 
equation near the horizon and then matching its solution to the solution near infinity. The spheroidal wavefunctions 

1) — s(s + 1). Matching the inner 



(|3.2p reduce to the spin- weighted spherical harmonics with eigenvalues s Ai m = 1(1- 
solution to the outer solution, the frequency ui is (see Appendix [A] for details) 



n 

k=l 



2iw 



2ivo 



1 - i(w(r + - r-)) 2l+1 L s (s + 2iw) IlLiO 2 + 4n7 2 - s 2 - Aiszu) 



1 + i(w(r + - r_)) 2i + 1 L s (s + 2iw) nLi( fc2 + 4ro2 



Aisw) 



(3.4) 



where 



- 1 (-i) s r(i + ;- s )r(i + ; + s ) 



T(l + l) 



T(2l + 2)T(2l + 1) 



(3.5) 



Equation p.4p can be solved numerically for the characteristic values of u>. Care must be exercised to ensure that the 
solutions are consistent with the approximation scheme, i.e., Mto <C 1, a <C M and zq ~ 0. The instability timescale 
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for scalar perturbations is 



iff = 2 



1 - (a/V2M) 2 + {a/Mf 



log" 



1 + 7(w„, m - mO) 



1 - 7(w„ 



(3.6) 



Tabic U shows the results for s = and different values of I = m. The scalar instability timescale is much smaller 

TABLE I: Approximate solution for a slowly rotating object with a = 0.3 M and mirror position at e = 10 -5 . The two columns 
give the real part of the frequency Re(uj) and the growth time r for different I — m values of the scalar field, respectively. The 
mode I = m = 1 is absent. 



a/M = 0.3, e = 10~ b 



/ = m 
2 
3 
4 



Re(uj)M t/M 

0.02 1.8 x 10 13 

0.09 5.8 x 10 13 

0.17 8.1 x 10 15 



than the Hawking evaporation timescale [261 ] . r e 



10 (M/Mq) sec, and increases with I = m. Gravitational 



instability timescales of fast rotating objects are expected to be larger by several orders of magnitude. 



B. Fast-rotating objects 



Defining x 
written as 



(r — r+)/r+, <7 = (r+ — r_)/r+, r = M(w — mO) and w' — » r+w, the Teukolsky equation can be 



[2iu'x 2 + x{AiJ - 2(s + 1)) - (s + l)cr + 4ir]i?'(a;) + [2(2s + l)iw'(a; + 1) + A]i?(.T) - a; (a; + a)B!'(x) = . (3.7) 

Equation (|3 . 7[) can be solved in the limit a —f M and u) — > mil by following the procedure of the previous section. 
The relation between the position of the mirror and the frequency of the wave is (see Appendix [B] for details) 



Ri + pR 3 {-2iua) 2t5 
R-2 + pR 4 (-2iujcr) 2i5 



(3.8) 



Numerical solutions of the above equation for a star with a = 0.998M are shown in Fig. [T] The real part of the 
characteristic frequency is always close to mfl and lm(w) <C Re(w). Thus the results are consistent with the initial 
assumptions. The instability timescale for gravitational perturbations is about five orders of magnitude smaller than 
the instability timescale for scalar perturbations. 



IV. INSTABILITY ANALYSIS: NUMERICAL RESULTS 



The oscillation frequencies of the modes can be found from the canonical form of Eq. (|3.ip 

d 2 Y 



dr 2 

where 



VY = 0, (4.1) 



y = A*/V+« 2 ) 1/2 ^ V - ^ ~ 2lS(r - ( V K + A(4 ^ S = A) - G 2 - f- , (4.2) 

and K = (r 2 + a 2 )uj — am, G = s(r — M)/(r 2 + a 2 ) + rA(r 2 + a 2 )~ 2 . The separation constant A is related to the 
eigenvalues of the angular equation by Eq. (|3.3[) . The eigenvalues s Ai m are expanded in power series of au> as [27j 



Terms up to order (aw) 2 are included in the calculation. Absence of ingoing waves at infinity implies 

Y - r~ s e lur * . (4.4) 
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Numerical results are obtained by integrating Eq. (|4.1[) inward from a large distance . The integration is performed 
with the Runge-Kutta method with fixed u) starting at Mr^ = 400, where the asymptotic behavior (|4.4[) is imposed. 
(Choosing a different initial point does not affect the final results.) The numerical integration is stopped at the radius 
of the mirror ro, where the value of the field Y(uj,rQ) is extracted. The integration is repeated for different values 
of u> until Y(u!,ro) = is obtained with the desired precision. If Y(ui,ro) vanishes, the field satisfies the boundary 
condition for perfect reflection and u) = loq is the oscillation frequency of the mode. 

A. Objects with a < M 

The regime a < M requires a surface or mirror at ro = r+(l + e) > r+. Typical results for scalar perturbations of 
objects with a < M are summarized in Table HI] and Fig. [2] The top panels of Fig. [2] show the imaginary and real 
parts of the fundamental mode frequency for a — 0.998A/ vs. the mirror position ro = r + (l + e) for different / = m 
values, respectively. The instability is weaker for larger m. This result holds also for I ^ m, although it is not shown 
in the plots. Comparisons with the analytic results in the near-extremal regime are shown in the bottom panels. The 
numerical integrations indicate that Re(w) ~ mfl, in agreement with the analytic results of Sect. IIII F3l The instability 
timescales are consistent with the analytic results within a factor ~ 3. The minimum instability timescale is of order 
r ~ 1 5 M for a wide range of mirror locations. 

TABLE II: Characteristic frequencies and instability timescales for a Kerr-like object with a = 0.998M. The mirror is located 
at e — 0.1. Results for scalar and gravitational perturbations are shown for different I — m values. 



(Re(w)M ,Im(Lu)M) 


I = m 


s = s = 2 


1 


(0.1120,0.6244 x 10" 5 ) 


2 


(0.4440 , 0.5373 x 10" 5 ) (0.4342 , 0.2900) 


3 


(0.7902,0.1928 x 10" 5 ) (0.7803,0.2977) 


4 


(1.1436,0.5927 x 10" 6 ) (1.1336,0.3035) 



Figure [3] shows the results for gravitational perturbations. Instability timescales are of the order of r ~ 2 — 6M. 
Thus gravitational perturbations lead to an instability about five orders of magnitude stronger than the instability 
due to scalar perturbations (see Tabic [TTJ) . Figure [3] shows that the ergoregion instability remains relevant even for 
values of the angular momentum as low as a = 0.6M. 

Some features of these results are intriguing and deserve further study. For instance, the regime Re(cj) > mf2 
may also be unstable. Since the superradiant instability is confined to the superradiant regime, this mode should 
describe a different kind of instability. Numerical results also show that there is a mirror location which maximizes 
the instability at fixed a. It would be interesting to explain the physical meaning of this location. 

B. Objects with a > M 

Objects with a > M could potentially describe superspinars. Several arguments suggest that objects rotating above 
the Kerr bound are unstable. Firstly, extremal Kerr BHs are marginally stable. Thus the addition of extra rotation 
should lead to instability. Secondly, fast-spinning objects usually take a pancake-like form [28| and are subject to 
the Gregory-Laflamme instability [29|, [3(J. Finally, Kerr-like geometries, like naked singularities, seem to be unstable 
against a certain class of gravitational perturbations [U HH, [33J called algebraically special perturbations dU . These 
perturbations are described by modes with zero Teukolsky-Starobinsky constant [23|], and will be discussed in more 
detail below. For objects with a > M the surface or mirror can be placed anywhere outside r = 0. In general the 
instability is as strong as in the a < M regime. An example in shown in Fig.[U which displays the resonant frequencies 
for the instability of a surface at tq/M = 0.001. This result seems to confirm other investigations suggesting that 
ultra-compact objects rotating above the Kerr bound are unstable (35j . 

C. Algebraically special modes 

Naked singularities are also characterized by a further kind of instability which is not related to the presence of an 
ergoregion [3l|, [32], [33| . This instability is generated by modes with vanishing Teukolsky-Starobinsky constant [23|. 
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FIG. 2: Details of the instability for scalar perturbations, from a numerical solution of Teukolsky equation. Top panels: 
Numerical results for the imaginary part (left panel) and real part (right panel) of the fundamental mode frequency vs. the 
mirror position ro = r+(l + e) for different / = m values. The angular momentum is a = 0.998A/. The instability grows 
monotonically with I = m. Bottom panels: Zoomed-in version of the upper panels, where numerical results (solid lines) are 
compared to the analytic solutions in the near-extremal regime (|3.8fl (dotted lines). The agreement between numerical and 
analytic values of the resonant frequency is remarkable. The analytic results for the imaginary part (bottom left panel) agree 
with the numerical results within a factor ~ 3. 



These "algebraically special modes" are regular for r > and may be relevant to superspinar geometries without 
pathological regions. They can be expressed in analytic form. The solution of the Teukolsky equation is 

Rlm(r) = (A + Br + Cr 2 + Dr 3 ) e 1 ^^-^"^ S lm (6) , (4.5) 

where A, B, C, D are constants and Si m are spin weight- 2 spheroidal harmonics [27j • These modes satisfy proper 
boundary conditions at infinity and are well behaved for any r > — oo. Although superspinars require particular 
boundary conditions at the excised region, unstable algebraically special modes can be present. They can be computed 
with the continued fraction method [3a. |37| and correspond to a zero of the Teukolsky-Starobinsky constant squared: 

|Cf = A 2 (A + 2) 2 - 8uj 2 \ [a 2 (6 + 5A) - 12a 2 ] + IUuj 2 (a 4 u 2 + M 2 ) , (4.6) 

where A is defined in Eq. (|3 . 3|) for s = — 2 and a 2 = a 2 — ma/uj. (Note a typographical error in Ref. [37|.) The 
technique of Ref. [36| can be used to evaluate the algebraically special modes at fixed a. In the range < a < M, the 
modes coincide with those of Ref. [37j • Some results for a > M are listed in Table IIIII The typical timescales are of 
the order of 10~ 6 seconds for a 1M@ star and 1 second for M = 10 6 M Q . 



V. DISCUSSION 



This paper presented a general method for investigating the ergoregion instability of ultra-compact, horizonless 
Kerr-like objects. The essential features of these objects have been captured by a simple model whose physical 
properties are largely independent from the dynamical details of the gravitational system. The method has been 
applied to superspinars and rotating wormholes. Numerical and analytic results show that the ergoregion instability 
of these objects is extremely strong for any value of their angular momentum, with timescales of order 10~ 5 seconds 




FIG. 3: Details of the instability for gravitational perturbations, for different I = m modes and a/M = 0.998 (top panels) and 
for I — m — 2 and different a/M < 1. 




1.14 1.16 1.18 1.20 1.22 1.24 1.26 1.28 1.30 1.32 




1.14 1.16 1.18 1.20 



1.22 1.24 1.26 1.28 
a/M 



FIG. 4: The fundamental I = m = 2, 3, 4 modes of an object spinning above the Kerr bound as function of rotation. The 
surface is located at ro/M — O.OOf. 



for a IMq star and 10 seconds for a M = 10 6 Af Q star. The above investigation confirms previous results for gravastars 
and boson stars (llj . namely that exotic objects without event horizon are likely to be ruled out as viable candidates 
for astrophysical hyper-compact objects. 
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2Mu 


M 


m=0 


m=2 


0.00 


(0, 4.0000) 


(0.0000, 4.0000) 


0.10 


(0, 4.03107) 


(0.918857,3.66121) 


0.20 


(0, 4.13218) 


(1.36523,3.07054) 


0.30 


(0, 4.33511) 


(1.51483,2.57005) 


0.40 


(0, 4.74699) 


(1.53973,2.19077) 


0.50 


(0.715878,6.57057) 


(1.51441,1.90331) 


0.60 


(2.47766,5.25615) 


(1.46928,1.68066) 


0.70 


(2.84578, 4.36141) 


(1.41732,1.50403) 


0.80 


(2.93170, 3.71616) 


(1.36428,1.36084) 


0.90 


(2.91142, 3.23046) 


(1.31270,1.24256) 


1.00 


(2.84632, 2.85259) 


(1.26369,1.14328) 


1.02 


(2.83051, 2.78693) 


(1.25424,1.12531) 


1.04 


(2.81414, 2.72409) 


(1.24490,1.10790) 


1.10 


(2.76243, 2.55079) 


(1.21763,1.05878) 


1.20 


(2.67192, 2.30455) 


(1.17458,0.98601) 


1.30 


(2.58073, 2.10007) 


(1.13444,0.92268) 


1.40 


(2.49185, 1.92767) 


(1.09701,0.86707) 



TABLE III: Algebraically special modes for various values of the angular momentum. 
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APPENDIX A: ANALYTIC SOLUTION IN THE LOW-FREQUENCY REGIME 

This appendix contains the analytic computation of the instability of rotating objects, bounded by a hard wall. 
Massless and massive scalar perturbations and general spin-s perturbations are considered. 

1. Massless scalar fields 

Following Ref. [2(1 HU, [25|, the space-time outside the star is divided in a near region, r — r + -C 1/w, and a far 
region, j — r+ 3> M. The radial equation (|3.1|) is solved separately in each of these two regions with the assumptions 
that the Compton wavelength of the scalar particle is much larger than the typical size of the object, l/u> 3> M, 
and a <C M. These solutions are then matched in the overlapping region, where the condition M <C r — r + <C l/u) 
is satisfied. The equation for the characteristic value w is obtained by imposing suitable boundary conditions at 
the boundaries. In the following, r + denotes the location of the "would-be" horizon. The location of the mirror or 
ultra-stiff wall is ro = r+(l + e), e <C 1. 

In the far region, where the effects induced by the BH can be neglected, one can approximate a ~ 0, M ~ and 
A ~ r 2 . The radial wave equation reduces to the wave equation for a massless scalar field of frequency ui and angular 
momentum I in flat background 

d 2 (rR) + [uj 2 - 1(1 + l)/r 2 } [rR) = . (Al) 

The most general solution of this equation is the linear combination of Bessel functions [38| 

R = r' 1 ' 2 [aJ i+1/2 W + /3J_ ; _ 1/2 (o;r)] . (A2) 

The behavior of Eq. (|A2[) for large r and small r are 

PT\ (iu/2) l+1 / 2 , (uj/2)-'- 1 / 2 , . 

R~sj-- [«sinK - l, t 2) +00^ + l«/2)] , R ~ a r> + , (A3) 
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respectively. Absence of ingoing waves at infinity implies (3 = — iae l7rl . The radial wave equation in the near region is 



Ad r (Ad r R) + r\(u - mft) 2 R - 1(1 + 1)A R = . 
Introducing the radial coordinate 

T — T-i- 

z = , < z < 1 , 



(A4) 



(A5) 



Eq. (|A4[) can be rewritten as 



z(\-z)d 2 z F + {(I + i 2m) - [1 + 2(1 + 1) + i 2m\ z) d z F - [(I + 1 f + i 2m(l + 1)] F = , 



where R = z itu (l - z) l+1 F and 



w — (uj — mQ)- 



r+ — r_ 

Equation (|A6[) is a standard hypergeometric equation. Its most general solution is 

R = Az-' L ™(\ - z) l+1 F(a -c+l,b-c+l,2-c,z) + Bz lzu {l- z) l+1 F{a, b, c, z) , 
where a = 1 + I + i 2m, b = I + 1 and c = 1 + i 2m. Near the wall, where r ~ r+, Eq. (j A8|) reads 

where r* is the tortoise coordinate 

r r 2 4- n 2 

dr . 



(r — r + )(r — r_) 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



The large-r behavior of the solution in the near region is obtained with the change of variable z — > 1 — z in the 
hypergeometric function (38j . The result is: 



1? 



r_i_ — r_ 
r 



r(2/ + 1) 
r(/ + i) 
-i-i 



„ T(l-2iw) T(l + 2im) 

A— A — - — —-- + B- 



T(l + l-2im) T(l + l + 2iw) 



r(-l - 21) 

n-i) 



, T(l-2iw) T(l + 2im) 

A —A- '- + B ■ 



Y(-l-2iw) T(-l + 2im) 



The matching of the near- and far-region solutions in the region M <C r — r + <C 1/lu yields 



B 
A 



R 1 +i(-l) l (m(r + -r_)) 2l+1 LR 3 
' R 2 + i(-l) l (uj(r + -r_)) 2l+1 LR 4 



where 



and 



L = 



_j[ (r(/ + i)r 

2 2/+ 2 F ( Z + 3/2)r(2/ + 2)r(2; + i)r(i + 1/2) 



r(i - 2iw) 



r(l + / - 2iro) 
Equation (I A12|) can be rewritten as 



r(l + 2iro) 
r(Z + 1 + 2im) 



r(i-2^) 



T(l + 2itu) 
r(-Z + 2iro) 



(All) 



(A12) 



(A13) 



(A14) 



1 



n( k + 2iro\ 
I fc - 2iro i 
fc=i v 7 



1 + 2L(r+ - r^) 2l+1 muo 2l+1 ]lLi( fc2 + 4 ^ 2 ) 
1 - 2i(r + - r_ f +1 vouj 2l + l ]lLi( fc2 + 4kj2 ) 



(A15) 
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where 

R% -y-r fk + 2im\ 

R~2~ ~ V/c-2in7 ) ' 



i? 3 = ii?i(-l) i+1 2tupj(fc 2 + 4zu 2 ), 



fc=i 

-ro have been used. If the mirror is located near the outer horizon at 



and similar relations for i?4 and i?2 with vj 
a radius r = ro, the scalar field must vanish at the mirror surface. This condition implies 



B 
A 



F(l + 1, 1 + 1 - 2zw, 1 - 2itn, zo) 
F(Z + 1, 1 + 1 + 2in7, 1 + 2izu, z ) 



(A16) 



where z$ = z(ro). The relation between the position of the mirror and the frequency of the scalar wave is obtained 
from Eq. (|A15|) : 



F(l + 1, 1 + 1 - 2iti7, 1 - 2iro, z ) 
F(l + 1,1 + 1 + 2im, 1 + 2im, z ) 



2ivj 



n 

k=l 



k + 2iw 
k — 2izu 



l + 2L(r, 



\2l+l 27+1 



nU(fc 2 +^ 2 ) 



1 - 2L(r+ - r-f l+1 wuj 2l + 1 IlLi( fc2 + 4 ^ 2 ) 



(A17) 



In general, Eq. (|A17[) must be solved numerically. However, an approximate solution can be obtained by assuming 
Re(w) 3> Im(w) and w -C 1, i.e. a frequency near the superradiant limit lo ~ mfl. This solution gives a good 
approximation for Mu> <§; 1, small m and slowly rotating objects. Since Re(w)/Im(oj) <C 1, Eq. (|A17[) can be first 
solved for real w, then a small imaginary part is added and the equation is solved again for Im(o>). The l.h.s. and the 
last two terms of the r.h.s. of Eq. (|A17[) are ~ 1 for frequencies near the superradiant limit. This yields zo 2m = 1. 
Using the tortoise coordinate, it follows 



a 2ivar°(r + — r_)/(2Mr + ) _ Hvjx 



1 



where r% — r*(zo) and x = r®(r+ — r_)/(2Afr + ) = log(zo). The solution of Eq. (IA18|) is 

7r(r + — r_) 



9 

r+x 



iVL. 



(A18) 



(A19) 



Positive frequencies can be obtained by imposing x < —nir(r + — r _) / '(mflri.) . The superradiant limit requires 
zu = nn(r + — r_)/(r?_a;) <C 1. This condition is satisfied by considering only the fundamental tone and the first 
overtones or placing the mirror very close to the horizon, \x\ ^> 1. By adding a small imaginary part to the resonant 
frequency, u> — w„ jm + iS, where S << u> n . m , Eq. (|A17[) becomes 



F(l + l,l + l~2i (ujq + ipS) ,l-2i(w + ipS) , z ) _ 2ivjQ 2pS 
F(l + l,l + l + 2i (w + ipS) ,l + 2i(m + ipS) , z ) ^ *° 



n 

k=i 



k + 2ivjQ 
k — 2ivoo 



1 + 2L(r+ 



\2Z+1 



^o<rnu(* a +w) 



1 - 2L(r + - r_ 



\2l + l 



(A20) 



where w = zuq + ipS and p = r3./(r+ — r_). The ratio of the hypergeometric functions in the l.h.s. of Eq. (|A20[) is 
~ 1 for the approximations used in the derivation. In terms of the tortoise coordinate, Eq. (|A20|) reads 



2pxS _ 



1 + 2L(r^ 



\2l + l 



^o<rnu(* a +w 



1 - 2L(r+ - r.f+^C lT fe =i( fc2 + W) 



21 + 1 Y\l 



(A21) 



The solution of Eq. (|A"2"Tj) is 



5 = Im(w) 



2r\x 



log 



1 + 7(cJ n , m - mfi) 



1 - 7(cJ„, m - mfi) 



(A22) 
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where 

7 S \r\ i^f Uf + to ° 2 ) r (i + 3/2) r(2l+wf +1 )r (i + 1 / 2 ) "" s ' < A23 > 

Both 5 and 0J n ,m are very small for ro ~ r+. However, the argument of the logarithm in Eq. (IA22|) is ~ 1 and the 
assumption Re(o;) 3> Im(w) is satisfied. 

The above results display two important features. First, Eq. (|A22|) and Eq. (|A23[) imply 5^0 for Re(w) ^ mfL The 
time dependence of the scalar field <E> is exp(— iuit) = exp(— iRe(w)t) exp(St). Thus the amplitude of the field grows 
exponentially and the resonant mode becomes instable for Re(w) < mf2. Second, x < x cr i t = — ri7r(?*+ — r-)/(mOr?_) 
and Re(tuo) oc nx~~ l imply 5^0 for rt 0. There is always a superradiant amplification for n > provided that the 
mirror position is closer than x cr u to the horizon and the approximations used in the above derivation are satisfied. 
The critical value x cl n is positive and outside the domain of x for n < 0. In this case the mirror can be located 
everywhere in the near region, but there is no superradiant amplification. The growth timescale for n > is given by 



t = 5 1 = 2px log 
or, in terms of the physical variables, by Eq. 



1 + 7(^n,m - mQ) 



1 - 7(w„ !m - mfi) 



(A24) 



2. Massive scalar field 

If the scalar field is massive, the wave equation is 

A Al A Al * = ^ 2 *, (A25) 
where the covariant derivative and fi is the field mass. The above equation is separable. The radial equation is 

d ( A dR\ + 2 2 + fl2 2 _ 4aMojr + +a 2 m 2 _ A , 2r 2 + ^Jl + ^ R = q ( A26) 

dr \ dr J 

Assuming au; <C 1, A ~ 1(1 + 1) and fir <C I, the solution in the near region is identical to the solution for the massless 
case. The equation in the far region is 

8 2 r {rR) + [to 2 -fi 2 - 1(1 + l)/r 2 ] (rR) = . (A27) 

The results of the massless case apply with the substitution u> — ► k, where k 2 = ui 2 — fi 2 > 0. The matching conditions 
are 



F(l + 1,1 + 1 - 2ivj, 1 - 2ivj, go) _ 2m J-r ( k + 2ivo \ 
F(l + !,/ + ! + 2im, 1 + 2iw, z ) ~ Z ° AJ- \k - 2iw ) 



1 + 2L(r + - tJ) 21+1 w(uj 2 - M 2 )' +1/2 Y[ l k=1 (k 2 + 
1 - 2L(r+ - r- f +1 m(u; 2 - M 2 )' +1/2 nU( fc2 + 4 ?n 2 ) 
The imaginary part of the frequency is identical to Eq. (|A22|) with 

21 i 



_ 2 / r + - T - 

7 = 7rr + 



Jlf + te ° 2 > Ui + y m I-!vm + m i + m } ( <" ~ p2) ' +V2 a ' (A28) 

The condition for superradiant amplification, ui < mfi, does not depend on the field mass. 

3. General spin-s fields 

Previous analytical and numerical calculations [25], have shown that superradiant effects for gravitational fields 
are stronger than for scalar fields. For instance, superradiant amplification factors are about 0.1%, 4.5% and 138% for 
scalar, electromagnetic and gravitational field, respectively. Since the effects induced by the BH in the far region can 
be neglected, the radial wave equation reduces to the wave equation of a massless field with spin- weight s, frequency 
ui and angular momentum I in flat background: 

rd 2 f + 2(1 + 1 - iwr)d r f - 2i(l + I - s)uf = , (A29) 
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where R — e lu>r r l s /(r). Introducing the radial coordinate x = 2iujr, the wave equation becomes a standard 
Kummer equation (38j. Its most general solution is a linear combination of confluent hypergeometric functions: 



R = e"^V _s (a M(l + 1 - s, 21 + 2, 2iur) +/3U(l + l- s, 21 + 2, 2iur)) . 
The large-r behavior of Eq. (|A30p is 



R 



r(2^ + 2) 
a r ( i + i— j E"") 



1-Z-s 



r l+2s 



a( _ 1)1+i - s r( 2 ; + 2) • 
1 j r(i + / + s ) p 



(A30) 



(A31) 



The first two terms in Eq. (|A31[) represent an outgoing wave at infinity and an incoming wave from infinity, respectively. 
The behavior for small r is 



R 



a + P 



r(-i - 2i) 



T(-l - s) 

Absence of ingoing waves at infinity implies 



J-s 



/3{2iuj) 



-1-21 



r(2/ + i) 



-l-l-s 



r(i + i-«) 



1 1 ' r(i + z + s ) 



The near-region behavior of the solution in the far region is 



R~a 



1 

— r" 
2 



(2i W )- i -^(-l) 



r(2/ + 2)r(2Z + i) 
r(i + /-s)r(i + z + s) 



„-l-l-s 



(A32) 



(A33) 



(A34) 



The radial wave equation in the near region is 



dr 



_ x dR 
dr 



K 2 - 2is(r - M)K 



A 



R = 0, 



where A = (I — s)(l + s + 1) + 0(luci). Using the approximate relations 

A'«r+-r_, K = uj(r 2 + a 2 ) - am w r\(uo - mfi) , K 2 — 2is(r — M)K w (r+ - r_) 2 



and introducing the radial coordinate z, Eq. (|A35p can be written as 

1-z 



z{\ - z)d 2 z R + ((a + 1) + (s - + 



+ (ro - is/2)' 



(l-s)(/ + s + l) 
1-z 



R = 0. 



(A35) 



Setting i? = z ltu (l — z ) l+s+1 F, the previous equation becomes a standard hypergeometric equation [38|]. Its most 
general solution is 

R = A z- s - lzu (l - z) t+1+s F(a - c + 1, b - c + 1, 2 - c, z) + B z 4CT (l - z)' +1+s F(a, 6, c, z) , (A36) 
where a = l + l + s + i 2m, b = I + 1 and c = 1 + s + i 2m. The behavior of Eq. (|A36|) near r + is 

R ~ J 4 e - s -'( w -" lSi )w r « + B e i ( w ~ mfl )OT r * . (A37) 

The large-r behavior is 



1? 



where 



r+ — r_ 



i?3 = 



r(2i + i) 



r(i + i) 

T(l - s - 2iw) 
Y{l-s + l- 2dm) 



(AR 1 +BR 2 ) 



r+ — r~ 



-l-l-s 



T(-l-2l) 



(Ai? 3 + Bi? 4 ). 



(A38) 



i?2 = 



T(l + s + 2im) 
T(l + s + l + 2im) 



i?4 



r(— l — s — 2im 
r(l + s + 2im) 



- 2 s 2 - Aism) 



rfl — s — 2im) , .. , -i , s -i-r , 9 

- i?i(-l)' +1 (.s + 2im) J[(k 2 + Am 2 

k=l 

I 

r , .,. , - J R. 2 (-l)' +1 (.s + 2 l ^)n(fc 2 +4 ro 2 - S 2 -4 isro ). 



(A39) 
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The matching of Eq. (|X38|) and Eq. (|X34|) yields 

-B _ _ tt / fc + 5 + 2iro \ 1 - ijLsjs + 2iw)g s fA4f)l 
A~ IJ^yk- s-2ivj) l + ijL s (s + 2ivj)g s ' 1 

where 7 = (iv(r + — r_)) 2 ' +1 , g s = Y\ l k= i{k 2 + izu 2 — s 2 — Aiszu) and L s is defined in Eq. (|3.5[) for / 7^ ±s — 1. Equation 
(|^40|) reduces to Eq. (|Xl~5|) for s = 0. 

If the mirror is close to the outer horizon, the field must vanish at the mirror surface. Using the small-r behavior 
of the solution in the near region and setting the radial field to zero at z = Zo, it follows A = —Bz^ +2z ^. The relation 
between the position of the mirror and the frequency of the wave is given by Eq. (|3.4p . 



APPENDIX B: ANALYTIC SOLUTION IN THE NEAR-EXTREMAL REGIME 

This appendix describes near-extremal horizonless objects by employing the analytic approximations of Refs. [25l . 
[39l | . The resonant frequencies for rapidly spinning BHs can be found by rewriting the Teukolsky equation in the 
Kerr incoming coordinate system {v,r,6,<p'). The latter is obtained from the Boyer-Lindquist coordinates with the 
coordinate transformation dv = dt+ (r 2 + a 2 )dr/ A, d(f>' = d<f> + adr/A. The radial equation for a field of spin s is (39j 

AR"(r) + 2[(s + l)(r - M) - iK]R'(r) - [2(2s + l)iur - A] R(r) = . (Bl) 

Defining x — (r — r+)/r+, a = (r + — r_)/?* + , r = M(u> — mfl) and u' — > t + uj, Eq. (|B1[) reads 

x(x + a)R"{x) - [2(2s + l)iu'{x + 1) + A] R(x) - [2iJx 2 + x{Aiuj' - 2{s + 1)) - (s + l)cr + +4ir] ii'(x) = . (B2) 

In the far region, x ^> max(er, r), Eq. (|B2|) reduces to 

x 2 R"{x) - [2iw'x 2 + x{HJ - 2(s + l))]R'(x) - [2{2s + l)iw'(x + 1) + X]R{x) = . 

Its most general solution is written in terms of confluent hypergeometric functions 



(B3) 



R(x) = ax- s -^ +2iu ' +i5 M Q + s + 2iJ +iS,l + 2iS, 2icux) + (3 {S 
where S 2 = Alu' 2 - (s + 1/2) 2 - A. The large-x behavior of Eq. (|B4|) is 



-S) , (B4) 



R~x~ 



aT(l + 2i5)(e™ /2iu}') 1 / 2+s+2iu} '+ i5 
T(l/2- s-2iu' + i6) 



„1— 4iw' 



aT(l + 2i5)(2iu r ) 



r\-l/2+s+2iw' -iS 



T(l/2 + s + 2iuj' + i8) 



+(3(6^-8). (B5) 



The term representing an incoming (outgoing) wave at infinity behaves as r 



— 1 — 2s ^liur* 



• /r for near-extremal BHs. 



(See Table 1 of Ref. [39(). Thus the first and second terms in Eq. (|B5j) describe an incoming wave and an outcoming 
wave, respectively. Absence of incoming waves at infinity implies 

T(l + 2iS) r(l/2 - s - 2vJ - iS) 2lS m _. 

Equation (|B4[) reduces to 

ii~ Qta .-«-V2+«"'+« + p x -B-lp+2bS-i6 ( B7 ) 

for small r. In the near region, Eq. (|B2j) reads 

x(x + &)R"{x) - [x{UJ - 2{s + 1)) - (s + l)cr + 4ir]i?'(x) + [2(2s + l)ia/ + A]i?(x) = . 

The radial wave equation in the near region becomes a standard hypergeometric equation with the substitution 
z = — x/er. Its most general solution is 



R = A z 1 - c F{a - c + l,b - c + 1,2 - c, z) + B F(a, b, c, z) 



(B8) 
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where a = 1/2 + s — 2iu>' + iS, b = 1/2 + s — 2iu>' —iS, c = 1 + s + in and k = —At /a. The mirror is located at r ~ r + . 
In this limit the solution behaves as R ~ Az 1 ~ c + B. The large-r behavior of the solution in the near region is 



R 



(- 



A(-l) 



t _ T(2-c)T(b-a) 



B 



T(c)T{b-a) 



r(6-c + l)r(l-a) T(b)T{c-a) 
xx-H A( _ 1}1 _ c T{2-c)T{a-b) 



or, in terms of the physical variables, 

-l/2-s+2itu'+i6 



R 



(;)' 



where 



i?i(<5) 



r(c)r(a - 6) 

r(o-c+i)r(i-6) 1 " r(a)r(c-6) 

r(2i*) [4 (-l) s+4K i?i + B R 2 ] + (S -> -5) , 

r(i — s — 



T(l/2 - s + 2iw' + iS)T(l/2 - 2iu>' + iS- in) ' 

r(l + s + in) 
r(l/2 + s- 2iuj' + i<5)r(l/2 + 2iu' + iS + in) ' 
R 3 = Ri(-5), R i = R 2 {-5). 



The matching in the overlapping region yields 



a = a b T(a - b) [A (-l) 1 " ^ + B R 2 ] , = a a r(6 - a) [A (-l) 1 " ^ + B i? 4 ] 



Using Eq. p6)). it follows 



where 



p{—2iu>'a) 



A(-l) s+lK R 1 + BR 2 



A(-l) s+iK R 3 + B R 4 ' 

T(-2i(5)r(l - 2i8)T(l/2 -s- 2iu' + id) 
T{2i5)T{\ + 2i5)T{l/2 -s- 2iuj> - iS) 



(B9) 



(BIO) 



(Bll) 
(B12) 

(B13) 
(B14) 



Equation (|B13[) reduces to Eq. (9) of Ref. for suitable boundary condition in presence of a horizon, i.e. A = 0. 
In the case under study, A is nonvanishing because the region containing the horizon is forbidden by the presence of 
the mirror. Thus Eq. flETg) yields 



B = +iK R x +pR 3 {-2iaLo') 2i& 
A 1 ' R 2 + pR 4 (-2i<ju;') 2lS 



(B15) 



Dirichlet boundary conditions at the wall require z s lh = —B/A. The relation between the position of the mirror 
and the frequency of the wave, Eq. IpTS]) . is obtained from Eq. (|B14| . 
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